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Abstract We study the symmetry classes for the random Dirac fermions in 2 
dimensions. We consider Nf species of fermions, coupled by different 
types of disorder. We analyse the renormalisation group flow at the 
order of one loop. At Nf large, the disorder distribution flows to an 
isotropic distribution and the effective action is a sigma model. 



1. Introduction 

Random Dirac fermions in 2 dimensions appear in various contexts in 
condensed matter physics. They were used to model disorder in systems 
with degenerate Fermi surface [2] , in the d-wave superconductors [3, 4] 
or to investigate the plateau transition in the integer quantum Hall effect 
[5] . In all these problems, the main question concerns the localization 
of quasiparticles, often different from the generic Anderson localization. 

The universal properties of a system where randomness is present de- 
pend largely on the discrete symmetries, like time reversal invariance 
or spin rotation symmetry. This principle governs the Wigner-Dyson 
classification of random hermitian matrices [1]. Considering other dis- 
crete symmetries, like particle-hole or chiral symmetry, led Altland and 
Zirnbauer [6] to extend this classification. However, it was recognized 
that the physics of 2d Dirac fermions is richer than expected from the 
random matrix classification, in the sense that for the same symmetry 
class, some systems may be critical and other not. One reason is that 
in two dimension, the sigma models, which are supposed to describe the 
universal physics, can support a topological term. In the form of a theta 
or WZW term, they can render the model critical [7, 8]. Sometimes, 
the presence of such a term can be deduced from general considerations 
[7, 8], in other situations it can be obtained exactly [3, 16]. A possible 
strategy, adopted in [14, 15] and applied here to all classes, is to try to 



1 



2 



derive the effective sigma model, topological term included, by chain of 
controlled approximations. The drawback is that we have to consider a 
large number of fermion species A'^^ in order to have a control parame- 
ter for the approximation. This means that the most interesting cases, 
where Nf is small, are out of reach of this method. 

The symmetry classes for 2d random Dirac hamiltonians, were re- 
cently analysed [9]. The constraint that the discrete symmetries should 
also preserve the form of the Dirac operator leads to a refinement of the 
Altland-Zirnbauer classification, in which the classes ALII, CI and Dili, 
appear each in two different forms. Here, we pursue the characterization 
of the symmetry classes for the Dirac Hamiltonians of [9]. We compute 
the beta function for all the disorder perturbations. When the number 
of fermion components Nf is large the flow of the coupling constants 
generically closes on a single coupling constant, a fact which simplifies 
considerably the derivation of the effective action. This was done for 
classes D in [14] and for classes C and CI in [15], and we refer to these 
papers for the details of the method. Without surprise, we obtain the 
sigma models defined on the manifolds identified by Zirnbauer [6] . The 
two different versions of classes ^III, CI and Dili give rise to sigma 
models with or without a WZW term. 

We consider random Dirac Hamiltonians in 2d, in the representation 

/ V + M -2id, + A-A 

\-2id, + A^ V-M ) ' ^ ' 

where dz = {dx — idy)/2, dg = {dx + idy)l2 and V, M, Az and Az are 
random matrix fields of dimension Nf. 

The discrete symmetries preserving the form of the Dirac Hamiltonian 
are [9] 

H = -PHP-\ P=(j^ PPt = l, p2 = l (2) 

H = ecCH^C-^, C=(^_['^ CC+ = 1, C^ = ±C.(3) 

The first condition will be called chiral symmetry, while the second one 
could be called particle-hole symmetry when ec = — 1 and time reversal 
symmetry (or spin rotation invariance) when ec = +1. The symmetry 
conditions impose constraints on the matrix fields V, M, A^ and A^, 
which are listed in [9] . We consider gaussian probability distribution for 
the disorder, compatible with the symmetry conditions. 

The spectral properties of the Hamiltonians 1 can be analysed using 
the supersymmetry [10] (or replica) method, in which products of Green 
functions are represented by gaussian integral over equal numbers of 
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fermionic and bosonic Dirac fields. Average over gaussian randomness 
introduce effective coupling between the bosonic and fermionic fields, 
the coupling constants being the disorder variances. 

The first step of our analysis is the renormalisation group flow for the 
effective coupling constants. The analysis can be eventually performed 
at all orders [11] . In the section 3 we give the one loop beta functions 
for all the clases and for any Nf. Generically, at large Nf the flow 
closes on a single coupling constant, corresponding to equal variances for 
different types of disorder. ^ In this case, the effective interaction, quartic 
in the the Dirac fields, can be decoupled via a Hubbard-Stratonovich 
transformation using a (super) matrix field Q whose dimension does not 
depend on Nf. The number of components of the Dirac fields appear 
as an overall factor and it can be used as control parameter for a saddle 
point approximation. 

The effective action for the field Q, hence the saddle point equation, 
has global invar iance, described by a group G. The solutions of the 
saddle point equation break the invariance to a subgroup H, so that the 
manifold of saddle points is given by G/H. The low energy behaviour 
is obtained by taking slowly varying configurations of the field on G/H 
(Goldstone modes) and evaluating the action for these modes. The result 
is a sigma model on some symmetric (super) space. Given the fact that 
in two dimensions sigma models can support a topological term, we may 
obtain a theta term or a WZW one. In the last section, we give a list 
of bilincars in Dirac fields to which couple the field Q, the invariance 
group G of Q, as well as a diagonal solution of the saddle point equation 
go together with its stabiliser H. Together with the coupling constants, 
these are the basic ingredients to characterize the sigma model for each 
symmetry class. 

2. Disorder perturbations 

We use notations similar to [14] . To obtain a compact notation, we 
introduce supermultiplets 



13+ 



</,* = (V'_, /?_, -01) 



(4) 



and similarly for , where ijjj^ and ^_ are fermion multiplets with Nf 
components while (3j^ and /3_ are bosonic ghostswith the same number 
of components. The superscript T denotes transposition. 
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In the absence of disorder, the action for the Dirac fields is written 

Scit = J ^{(l>%ct^ + ^%^), (5) 

and it corresponds to a conformal field theory with central charge c = 0. 
Prom the holomorphic (antiholomorphic) fields (f> {4>) > one can construct 
the current algebra OSp{2Nf\2Nf)k=i [16] . 
The perturbation induced by disorder is 

Spert = 9iJ^Oi, (6) 

where gi are the disorder variances associated to different disorder in- 
variants and the operators Oi are defined as follows 

leA _ 
lev 

Om = ^(<^%+</')^ Om = {^'E+(t>)\ (7) 
I&M 



where 



„± _ fei 



ej are generators of gl{Nf) normahzed by Tr(e/ej) = 6ij and A, V, M. 
are subsets of indices of sl{Nf) which depend on the symmetry class, eo 
is the identity normalized by TV Cq = 1. 

3. The one loop beta functions 

To study the renormalisation group flow for the coupling constants, 
we need the corresponding beta function. They can be derived from the 
operator product expansion of the perturbing operators Oi 

Oi{z, z)Oj{0) ~ ^Cf,Ofc(0) + reg. (8) 

The one loop beta functions are given by [13] 

Pk = ldigk = -Y.ctjgigj (9) 

The results for different classes are listed in the remaining of this section. 
The non-chiral classes 
Class {A) 
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This class do not posses neither chiral invariance, nor particle-hole/ 
time reversal symmetry. We set Nf = N and consider only disorder which 
do not break parity invariance of the action. The disorder fields A, 

V and M belong to gl{N); the sl{N) part generates the perturbation 
operators Oa, Oy, Om, while the gl{l) part generates Oa, Oy, Om- 
The corresponding beta functions are 

D 2 D f ^ 2\ B 1 1 

PA = --^9 A - [gu + gv) + -^gMgv + -j^gmgv + j^gngv , 

Ny 1 

Pa = —^gugv + j^gmgv , 

Pm = -xgii - ^ (5m + gv)gM + {^x - ^ gAgu + ^gAgv 

2 2 

+ j^gAgv + j^gaigv + gn) - xgugv , 

2 1 / ^ f D\ B 

pv = xgy + —{gm + gv)gv + I ~ y ) 9Agv + -^gAgu 

2 2 

+ j^gAgm + j^gaigM + gv) + xgngv , 

= -j^gm{gm + gv) - CgmigM + gv) + '^CgAgm 

+ j^gaigm + gv) + -^dAgv , 

Pv = j^gv (gv + gm) + Cg^ {gv + gu) + "^CgAgv 

2 . , Ny 

+ j^gaigv + g-m) + -^gAgn , 

where the constants, with group theoretical significance, are 

_ - 1 2(N^ - 4) 

1 _ 4(7V2 - 1) 

In the large N limit, the line = gy = gu = 9,n = -Qv = "2^9 a = 
g ~ \/N is preserved by the RG fiow. It is interesting to note that we 
cannot obtain an invariant line with all the coupling constants positive. 

Classes 3_ and 4_ {D and C) 

These classes have particle-hole symmetry {e^ = —1), inforced by 
symmetric, respectively antisymmetric matrix = zta. We set = 
2N, but for class D, N can be half integer. The random fields A and 

V belong to the algebra so{2N) (case D) and sp{2N) (class C). These 
algebra are spanned by a subset of the sl{2N) generators e/, r/, selected 
by the condition tj = —aTja~^. The field M belongs to the complement 
of these algebras with respect to sl{2N), spanned by Tj with Tj = 
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aTja ^. Operators Oa, Oy are not allowed by these symmetry classes. 
The beta functions are 

PA = — — -^9v + ^^9Mgv + -^9m9v;, 

Pm = -XtTT9m - -^9m9M + {^XrTr - "^^^ 9AgM + ^^gAgv 
— XrTrgMgV , 

PV = Xrrrgv + 1^9m9V + y^Xrrr \ gAgv + —^gAgu 

1 

+ j^gmgA + XTrTgugv , 

Pm = —^9m + '^CgA9m + Ny^rgAgv " Cgygm - C'gmgM , 

with the group coefficients (again, the upper sign is for class D and the 
lower one for class C) : 



C 



2 — ' - ^sl{2N) - ' 

Drr = -2(7V:f1), D^t = -'^N , Dtt = -2{N ±1) , 

1 l^N l±N 
Xttt = —XtTt = ='=2 ' XTTT = , XTtT = ^jy 

Vrr = Brr = 2iN^l), 

R o ^'-l R (Ar±2)(iV^l) 



N ' N 

At large N, an invariant line is given by gA = gM = gv = gm = g- 
Classes 3+ and 4+ (All and AI) 

These classes have time reversal (spin rotation) symmetry, with ec = 1 
and a'^ = iter. They are similar to classes 3_ and 4_, with the random 
fields V and M changing roles: now M belongs to so{2N) and sp{2N) 
respectively, while V belongs to the complement with respect to sl{2N). 
The beta functions can be obtained from that of the classes 3_ and 4_ 
by changing gM •s-^ -gv and gm ^ -gv 

The chiral classes 

• The first type of chiral classes corresponds to 7 = 1 in eq. 2 and 
it comprises three classes. Only disorder of the vector potential type 
is allowed, and these classes can be seen as limits of the previous ones 
when^y, gM, gv, gm ^0. 

Class 1 {AIII) 

The random gauge potential A belongs to gl{N) 2± gl{l) sl{N) 
{Nf = AT). 

PA = -jgl, Pa = 0, (10) 
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with D = -2N. 

Classes 5 and 6 {Dili and CI) 

In addition to the chiral symmetry, these classes posses particle-hole 
symmetry. We set Nf = 2N. The random gauge potential A belongs to 
so{2N) and sp{2N) respectively 

Pa = -^9\ , (11) 

with Drr = -2(iVT 1)- 

• The second type of chiral classes correspond to 7 = o"3 and it com- 
prises five classes. We set Nf = 2N, but for classes 8 and 9_|_ N should 
be understood as even integer. 

Class 2 (AIII) 

The chiral symmetry imposes the constraint [A, 7] = to the ran- 
dom gauge potential; this condition selects an algebra gl{N) ® gl{N) ~ 
{gl{l) (g) sl{N))®'^ out of gl{2N). The two gl{l) parts are generated by 
the identity and 7 and we denote the associated coupling constants by 
ga and ga-- We denote by ej the generators of sl{N) (g) sl{N). The 
two sl{N) components can have different variances, gAi and gA2j with 
5a = i9Ai ± gA2)l2. 

The fields V and M are constrained by {1^,7} = 0, {M, 7} = 0, so 
they belong to gl{2N) \ {gl{N) (g) gl{N)), generated by e^. Operators Oy 
and Ojn are not allowed by the chiral classes. 

The beta functions are 

Pa = — ^ial +9a) ^ [9m +9v)+ -^9Mgv , 



a- 7~i + - /3 ^Ua n Daa / 2 , 2 A 

Pa = -DugAgA ' Po. = -^9Mgv , Pa- = — ^ [gn +9v) ^ 

Dirv _L Birv _L 1 



Pm = -Xaaaguigu + gv) — y3a9m + -^gA9v + j^gaigM + gv) , 

D' B' 1 

Pv = Xaaagvigv + gu) — ^gXav + -^gtgM + j;^ga{gv + gu) ■ 

with the group coefficients 

Da = -Ds/(Ar) = —2N , Daa = Dgi{2N) - -DpZ(Ar) = -2iV , 

Dia = Dgl{2N)/'^ = -2iV + -| , = 2 , 

- R _2(7V2-1) o _2/v 

At large iV, the line g~X = 0, g\ = gM = gv = '^ga = ^ga- is left 
invariant by the RG fiow. 

Classes 7 and 8 {Dili and CI) 
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These classes have both chiral symmetry and particle-hole symmetry, 
with the matrices inforcing these symmetries commuting. 

Wc split the generators of gl{2N) in four groups Tj, Tq,, Tj, Tq, with 
Tj, Tj commuting with and Ta, anticommuting with it; also, 
"^i.a = ~^'^I,a^~^ aiid Ti^a = (^T^aCF^^ . The random gauge potential 
A is an element of so{N) ® so{N) or sp{N) fS> sp(A^)generated by Tj. If 
we choose the representation €c = —1, V belongs to o{2N) \ o{N) i^o{N) 
or sp{2N) \ sp{N) ® sp{N) respectively, generated by Tq, while the field 
M is generated by Tq. 

The beta functions are 

(3-A = -DiUhJ, 

T)ii T^OiOi 7~)aa oaa 

Pa = ^IfA +9a ) ^9m ^9v + -^9Mgv , 

(jjia \ nia 
^^tTt ^ j 9a9m + —^3a9v - x^TT9Mgv , 

^^TTT ^ I S^^y + —^9a9M + XTrT9Mgv , 

where the group coefficients are 

Di\ = -N±2, D^^ = D^'^ = -N , 
D\% = Di^ = -N±l, 

aaa _ aaa _ ^aaa _ aaa i ■'^ iai _ iai _ n 

B^^ = N, s;^ = B;^ = Ar:f 1. 

In the large N limit the RG flow preserve the line = 0, g'^ = gu = gv- 
Classes 9+ and 9 (Dl and CII) 

These classes have also both particle-hole and chiral symmetry, with 
the matrices enforcing these constraints anticommuting, {7, a} = 0. 

If we choose ec = — 1 , the matrix a defining the particle-holc symmetry 
is symmetric for class Dl and antisymmetric for class CII, as for the cases 
3_ and 4_ or 7 and 8. The opposite choice exchanges the symmetry of 
the matrix a between the two classes. 

As in the previous case, wc split the gl{2N) generators into four 
classes. Since {o'l'^} — 0) '^i spans now gl{N), to which the field A be- 
long. V is generated by Tq,, which span o{2N) \gl{N) and sp{2N) \gl{N) 
respectively. M is again generated by Tq. 

da- = -^g'u-^g'v, P2 = -D%gh-A, 

T~)ii r)aa T-taa iDaa 

Pa - — 7ry9A +9a) — i— fM — ^gv + —r9Mgv , 
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PM = -xttt9m + I -^x^tt Y' } ^^^^ ~^9a9v - x^TTdMgv 

Pv - x^^^ gy + 2x^^^ — \ Qj^gv + —^gAgM + Xj^^^gMgy , 



with 



^rT - ^ , ^rT- 

aaa aaa aaa aaa n ^iai ^ ^ 

Irprpj, — X^^^ — X^rp^ — J^^tT " ^ ' "^tTt ~ 27V" 2 ' 

1 1 

— ± - 

2N 2 



At large A'^ the RG flow preserves the line 5^ = and g^ = gu = 
gv = ^ga-- 

We note that the equations listed above have interesting properties 
under the transformation qm ^ —gv and gm ^ —gv This operation 
leave the equations for classes 0, 2, 7 and 8 invariant and exchanges the 
ones for classes 3±, A± and 9±, therefore it is equivalent to changing 
the sign of ec, when defined ^. Moreover, on the line qm = —gv = g 
(= g-m = —gv, when they exist), the flow for g and the flow for qa 
decouple. This was noted in [18] for class and was related to "spin- 
charge" separation, for class C in [12] . 

4. Hubbard-Stratonovich transformation and 
the sigma model 

We consider products of n Green functions, ur retarded and ua ad- 
vanced, and we need to introduce a copy of the multiplcts 4 for each 
Green function. When the Hamiltonian has chiral or particle-hole in- 
variance, the spectrum is symmetric with respect to £^ = and we do 
not need to distinguish between retarded and advanced sector. In this 
case, the density of states can be singular at the point E = Q. 

As already mentioned, when the coupling constant have equal values, 
the decoupling of the quartic interaction is particularly simple. At that 
particular point, the disorder perturbation can be written only in terms 
of singlets of gl{Nf). The decoupling matrix Q couples to these singlets 
and it has a size independent of Nf. Up to an additive constant, the 
effective action for Q is 



S[Q] = -Nf 



I ( ^^STVQ2_ STrlnf 
g J 2iT V 



Q d 
B Q 



(12) 
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where g ^ Njg. Part of the invariance of the free action 5 , OSp(2nNf\2nNf), 
survives in the effective action 12. We denote the residual invariance 
group by G. 

The saddle point equation 5S[Q]/5Q = can be solved making a 
diagonal ansaz Qo = /J- qo, where = 1 and /x has the dimension of 
a mass. As there is no scale in the problem, gi being dimcnsionless, fi 
is dynamically generated. The group G acts on Q by conjugation, and 
from a particular saddle point solution qq one generates a whole manifold 
by 9o q = gqog~^- This manifold is isomorphic to G/H, where H 
is the stabilizer of go- The last step in deriving the sigma model is to 
evaluate the action 12 on slowly varying configurations q{x) on G/H. 
An elegant way of doing it is to use non-abelian bosonisation [17], as in 
[4, 14] and to write the free action as a WZW model . The field q{x) 
becomes a mass term for the WZW field, which is forced to follow it. At 
scales larger than the eff'ective action is 

S[q\ =-Y^J^ STr d,qd,q + StoM (13) 

In the following, we give the main characteristics of the sigma model 
manifold for the 13 classes. 
Class {A 

Q ~Trgi(N) [S3,(/>,^* + ,^(/.*] , go = S3® A, 

G = GL{n\n) , H = GL{nR\nR) (g) GL{nA\nA) , 

with A = diag(l„^, -In^). 

Classes 3_ and 4_ {D and C) 

Q ~ IVgi(2N) S + , go = S3, 

G = OSp{2n\2n) , H = GL{n\n) , 

with E = diag(l, a) 

Classes 3+ and 4+ (All and AT) 

Q ~ IVgi(2N) S (S3#* + #*S3) go = A, 

G = OSp{2n\2n) , H = OSp{nR\nR) (g) OSp{nA\nA) ■ 

• The effective action for the classes 1, 5 and 6 is invariant under 

multiplication at left /right by holomorphic/antiholomorphic group ele- 
ments. This is a sign that the theory is conformally invariant and indeed, 
when deriving the sigma model we obtain a WZW term. 
Class 1 (AlU) 

Q ~ TVgi(N) [S3, #*] , 0-1/2 IVgi(N) [E3, #*] , 
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qo =1:3(8) A, G = GL{n\n)L (S: GL{n\n)R , H = GL{n\n) . 

The level of the WZW action is k = Nf. 
Classes 5 and 6 (Dili and CI) 

Q ~ Trgi(2N) (S</)<^*S-1) , Q ~ Trgi(2N) (s0</)*S-i) , 
qq = S3 , G = OSp{n\n)L 'S' OSp{n\n)R , H = OSp{n\n) . 

The level of the WZW action is fc = zb2A^, the change of sign indicating 
that the compact and non-compact sector switch from the orthogonal to 
the symplectic sector from one class to another. 

• To take into account the chiral symmetry for the last five classes, 
we double the supermultiplets by 



V2 \lt) ' V2 

Class 2 (AIII) 



Q --Trgi(2N) P3,^^* + , go = S3®f^3, 

G =GL{n\n)®GL{n\n) , H = GL{n\n) , 

where o"3 is the Pauli matrix acting in the extra space introduced to take 
into account the chiral symmetry. G is selected from GL(2n\2n) by the 
condition [5,cri] =0. 

Classes 7 and 8 (Dili and CI) 

Q ~ T>gi(2N) S + , go = S3 (8 (73 , 

G = OSp{2n\2n) (8 OSp(2n|2n) , H = OSp{2n\2n) . 

G is selected from OS'p(4n|4n) by the condition [g, ai] = 
Classes 9_ and 9+ (DI and CII) 

Q ~ IVg^sN) S + S-^ go = S3 (8 (73, 
G = GL(2n|2n) , H = OSp{2n\2n) . 

The group G is selected by from OSp{An\4n) by [^, S3 (8 (7i]. 

The classes 2, 7, 8 and 9-i- do not allow a theta term; also, there is no 
WZW term. Note that the sigma model is defined on the same manifold 
for the pair of classes 1 and 2 (class ^III), 5 and 7 (class Dili) and 6 and 
8 (class CI) ; the difference between the two members of a pair consists 
only on the presence or absence of the WZW term. 

The coupling constant for the sigma model is given by the inverse of 
the number of the Dirac fermions, / = ^/Nf. 
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Notes 

1. Exceptions may be the coupling constants corresponding to some u{l) part of the 
disorder, specific to classes 0, 1, 2 and 9±. This situation was analysed in the case 1 (AIII) 
in [4] and in [18] in a non-hermitian version of case 0. For the moment, we ignore this point 
and analyse the case of equal disorder variances. 

2. We thank to D. Bernard for pointing out that this can be realized on the Dirac Hamil- 
tonian via a non-unitary transformation. 
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